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The quantum noise in a linear amplifier is shown to be thermal noise. The theory of linear 
amplifiers is applied first to the simplest, single or double oscillator model of an amplifier, and then 
to linear model of an amplifier with continuous fields and input and outputs. Finally it is shown 
that the thermal noise emitted by black holes first demonstrated by Hawking, and of dumb holes 
(sonic and other analogs to black holes), arises from the same analysis as for linear amplifiers. The 
amplifier noise of black holes acting as amplifiers on the quantum fields living in the spacetime 
surrounding the black hole is the radiation discovered by Hawking. For any amplifier, that quantum 
noise is completely characterized by the attributes of the system regarded as a classical amplifier, 
and arises out of those classical amplification factors and the commutation relations of quantum 
mechanics. 



I. INTRODUCTION 



5h . 

Linear amplifiers, devices which take in a signal and produce and output signal of a different amplitude, are 

ubiquitous, but in general seem to be poorly understood. All produce noise, but again the source of that noise tends 

to be poorly understood, and seems often based on a case by case analysis. While often an amplifier can have excess 

noise, caused by some infelicity in its construction, all amplifiers must have a minimum level of noise, set by quantum 

Q\ ' mechanics. This was recognized by Haus and Mullen[l| and other 0H 00 na l r a century ago, but the lesson bears 

^sO , repeating. In particular, the noise is often roughly characterized by a temperature. We will see that this is exact- the 

noise output is thermal. Furthermore, black holes turn out to simply be an unusual instance of one of these amplifiers. 

By a linear amplifier I mean a device into which one feeds an input signal I(t) and out of which comes an output 

£ — . ' signal J A(t — t')I(t')dt' . In principle, that amplification A could be a function of both t and t' rather than just their 

difference. Such an amplifier is in general a phase sensitive amplifier, for example one which amplifies the cosine 

^ . component of the input differently from the output, but I will here be interested in phase insensitive amplifiers as 

above. 

As the simplest model of an amplifier let me first assume that the input and the output are both single modes, 

defined by the input canonical variables x, p and output Y, P. They are not time dependent signals- they are simply 

single degrees of freedom, which change in time only due to the dynamical Hamiltonian for that degree of freedom. 

C^ ■ The only "signal" is the value of that variable. An example of an amplification would be that the output degree of 

freedom, represented by Y and P is A times the input. 

Y = Ax (1) 

P = Ap (2) 

No matter what the input, the output is A times larger. However, it is clear that this amplifier is unrealizable. Any 
amplifier is some physical device which produces a unitary transformation between the input and the output. In 
particular, if Y and P are conjugate variables, we have 

[Y,P] =A 2 [x,p] = A 2 ih (3) 

But clearly this is the correct commutation relation only HA 2 — 1, and the amplification is trivial. In order to have 
such an amplifier one cannot have only one input channel. One needs at least two. (In this case a "channel" just 
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means another degree of freedom.) 

Y = Ax + Bq (4) 

P = Ap + E9 (5) 

where is the conjugate momentum to q. Demanding that [Y, P] = ifi then leads to 

A 2 + BE = 1 

We can always do a canonical transformation of the form e ri q and e~ v 9 to have 

E = -B (6) 

so that 

A 2 - B 2 = 1 (7) 

or A — cosh(/i) and B = sinh(/i). 

Clearly if one has two input channels, one also needs two output channels as well. Let me designate the second 
output channel by Z, P. Then in order that the commutation relations of Z and P and Y and P be maintained, we 

need 

Z = Aq + Bx (8) 

P = AO - Bp (9) 

(One could also have other canonical transformations of Z, P which would of course leave the commutation relations 
the same, but this is the simplest case, and using one of the others do not change the results). 
This is more easily expressed in terms of creation and annihilation operators Defining 



ip 



V2 ' 



(10) 



we have 



° - *# 

C = cosh(fi)a + sinh(/i)6 t (14) 

D = cosh(^)6 + sinh(^)a^ (15) 

or 

a = cosh(n)C - sinh( J u)D t (16) 

b = cosh(/i)L> - sinh(/x)C t (17) 

This is the form of a Bogoliubov transformation. 

Note that the commutation relations are maintained if we multiply a or b by phase 

C = cosh(iA)e iv a + smh(n)e- iK tf (18) 

D = cosh(^)e lK b + smh(fi)e- w a^ (19) 

or multiply each term overall by a phase factor. In the following I will not follow this complication, since all it does 
is to make the equations messier. 

Let us now assume that the input states of the two input modes, represented by o, b are thermal states, with input 
density matrices 

Pa = N a e- A « afa (20) 

Pb = N b e A ^ b (21) 



N a and N b are normalisation factors equal to 

N a = 1 - e- A » (22) 

and, if a has frequency w a , then A a = ^ where T a is the temperature of the a input channel. 
The density matrix in terms of the output annihilation and creation operators C, D is then 

p = N a N b e- A * afa e- AbbJib (23) 

= N a N b e' A ^ a ' AbbH (24) 

at at — A a (cosh(^)C^ — sinh(/i)D)(cosh(^)C— sinh(fi)D^)— A b cosh([i)D^ — sinh(fi)C)(cosh(fi)D— sinh(/^)C^) (OK} 

= N a N b e~( Aa ^ cosh2 (M)-A(, sinh 2 ( M )(C t C) + (A a +A b ) cosh( M ) sinh( M )(CD+Cf,Dt) (26) 

Taking the trace over the D using the complete set of m particle states D^D\m) = m\m) 1 X)m( TO H TO ) an d 
expanding the exponential in a power series, we see that each term must have the same number of D as D^ operators. 
((m\D' r D s \m) ~ Ounless r = s). Since each D in the expansion is multiplied by cither a fl' or a C and D^ by 
D or Cf , one must thus also have the same number of C as Cf operators. The reduced density matrix is thus a 
function only of C^C (after an appropriate number of commutations). It is in fact an exponential function in C^C 
(see appendix) of the form e~ ° c , where we can find Ac from 

TrC^CN c e~ AcC " c = -N c d Ac Tr(e- AcC " c ) (27) 

p-Ac 

Thus 

-Ac 

- Tr(iV a 7V 6 cosh(^) 2 a t asinh(/x) 2 (b t fe+l)+cosh(/x)sinh(/x)(a6 + a t 6 t )e- A « ata - A66t '') (29) 



1 - e- A c 

\2 



. £^ + sinh( M ) 2 ( * +1) (30) 

since the a, b modes are uncorrelated. For small A ab we get 

1 cosh(/i) 2 , . .,. 1 , . 

T~= A +sinh(^) 2 (— + 1) (31) 

Ac A a Aft 

If the inputs and outputs all have the same frequency io then A = ^ and this becomes 

T c = cosh(fi) 2 T a + smh(fi) 2 (T h + to) (32) 

and if the amplification is large, so that cosh(^i) w sinh(/i), 

-L=cosh( M ) 2 (-^ + -L) (33) 

Ac A a A b 

or 

Tc = cosh( M ) 2 (T a + T b ) (34) 



(recall that cosh(/i) is the amplification factor A). 

For A aj (, large (which corresponds to low temperatures) , we have 



e A ^ -1 
or 



= cosh(p) 2 e~ Aa + sinh(^) 2 w sinh(^) 2 (35) 



A c = -2in(tanh(^)) (36) 



This is exact in the limit as the input temperatures go to zero (A a- f, — > oo). Ic, for low temperatures in the inputs 
(temperatures much less than the input frequencies) , the output temperature is determined by the amplification and 
the frequency of the output solely. This is a quantum noise. Assuming the output has frequency ujc we have 

T<J = -2/n(tanh(//)) ^ 



Note that the output thermal noise due to the amplification of the vacuum fluctuation is given purely by the 
amplification cosh(/i) and the frequency of the output, both of which are determined purely by the classical behaviour 
of the amplifier. The quantum behaviour of the amplifier is a "classical" effect, in that it depends only on the classical 
attributes of the amplifier. That the expression for the temperature includes both r- does not alter the fact that it 
is completely determined by classical measurements. 

If we put the input into a coherent state, a\a) — a\a), b\a) = 0, then 

(a\C\a) = cosh(^)a (38) 

where a can be as large as desired. Ie, by measuring the output for a classical input, one can determine the parameter 
of the amplifier which determines the noise output of the amplifier. 

Alternatively one could have a situation in which one takes D as the output channel to be measured with a still 
being the input channel. Then the amplification of a in the D output is 

(a\D\a) = sinh(X)a (39) 

I.e., for small /i the "amplification" goes to zero, rather than to 1. 

It is also of interest to note that while the two inputs are, by assumption, statistically independent (no correlations 
between a, b) the outputs are not. Even in the vacuum input, we have 

(CD) = (0| cosh(^) sinh(ju) (aa f + 6+6)10) = cosh(^) sinh(^) (40) 

which implies a correlation (entanglement) between the Cand D outputs. 

That same entanglement implies that we could have "noiseless" (ie, not altering the signal to noise ratio of the 
input signal) by choosing an input state which was an entangled state- ie such that in the output state, the C and D 
modes were in a product state[6|]. One would then have a noiseless (zero temperature) output. This is in general not 
possible, because the input signal comes from somewhere that one is not able to correlate its quantum fluctuations 
with the other input to the amplifier. However in certain situations, in which the signal is a classical signal imposed 
on a quantum input channel, this may allow one to reduce the noise in a detector by choosing an appropriately 
correlated set of input channels, as for example in and interferometric gravity wave detection in which the gravity 
wave signal, a very large "classical" source, affects a quantum input channel in the electromagnetic field in the arms 
of an interferometer. 

II. CONTINUUM 

While the above "two mode" analysis is important, it is also instructive to examine a model for a continuous phase 
insensitive amplifier- ie, one with a continuous input and an continuous amplified output. 
Let us define the Lagrangian 



L= l - 
2 



2 - {d x <pf - (V? - (d^) 2 ) + 2q(# + i^)6(x - A)J dx (41) 

+ \tf+<»V) (42) 

with reflection boundary condition at x — of d x ip(t, 0) = d x (f>(t, 0) = 0. A is assumed to be very small, and we will 
take the limit as A — » 0. 

I could have taken x to be a continuous variable with the field propagating from infty to +oo and the oscillator 
located at x — 0, but in that case all of the antisymmetric modes for the <j> and ip fields would not have interacted 
at all with the oscillator. I take the oscillator at A and take the limit as A goes to zero to ensure that I am properly 
treating the interaction. The 6 at x = is not well defined on the half line x > 0. 

The ip field has negative definite energy, which is the source for the energy amplification which accompanies the 
amplifier. (Note that while this particular model amplifies energy that is not necessary for an amplifier, as we will 
see below) 

The equations of motion for the field are 

d 2 t 4>-d 2 x <t> = ed t qS(x-X) (43) 

dfy-d^ ~- -edtqS(x-X) (44) 

}2„ , r.2 



dfq + Wq = -ed t cf>(t,\)-ed t tp(t,X) (45) 



which have solutions 






e 2 - ? 2 



9 2 g + Q 2 g + — ^ — at(ff(*) + q(t - 2A)) = -ed t (Mt - A) + 0o(*+ A)) - e$(^(t - A) + ^ (* + A)) (48) 

where I will only be interested in the limit as A — > 0. Taking the Fourier transform of the resulting equations where 
q(t) = J q^e^ 1 ^, we have 

. (#oM + gW) ^o\ 

&> = 2zw 2 ~~n. ~2\ i r>2 (49) 

— or — iu){e z — e J ) + S2 J 

We take <po{t + re) and tpo(t + x) as the ingoing modes (the x and q of the above simple two mode analysis) and 

(pout = ipo{t-x) + eq(t-x) (50) 

tpout = ipo(t-x) - eq(t-x) (51) 

are the outgoing modes corresponding to C, D of the simple two mode analysis. 
Thus we find 



e4> (uj) + eip (uj) 
-oj 2 - i(e 2 - e 2 )oj + D, 2 



o „t(u;) = Vo(w) + 2»^e ■ _ _ (52) 



The conserved norm for the system is 

f 9 t - i>'*d t il>)dx + q'*(d t q + e<f>(t, 0) + c^(t, 0)) (53) 

d t O - d t i>'*i>)dx + {d t q'* + e^*(t, 0) + e^*(t,0))g 



(54) 

where S = {cj>, -0, 5} designates a complete solution of the equations of motion at any time t. This norm is conserved 
by the equations of motion and relates the ingoing modes at t — > — 00 to the outgoing at t — ¥ +00. 

Note the sign of the ip term in the norm. This arises from the fact that the conjugate momentum for the tp field is 
-dti>. 

The quantization of the the fields is such that positive norm fields are associated with annihilation operators while 
the negative norm fields are associated with creation operators. In the case of the ip field, the vacuum state, annihilated 
by the annihilation operators 

a|0) = (55) 

is also a maximum energy, rather than a minimum energy state. Also, while as usual, the positive norm states for 
the 4> fields are the positive frequency states, e~ lu)t the positive norm states for the ip field are negative frequency 
states e lult . Thus, the outgoing positive norm <f> states are linear combinations of the ingoing positive norm (j) states, 
and ingoing negative norm ip states, and the annihilation operators of the outgoing <p> field are linear combinations 
of the annihilation of the ingoing <j> field and creation operators of the ingoing ip field. This is precisely the situation 
examined in the first section. 

The fact that the ip field has negative energy is clearly an approximation in any real world situation, as the energy 
will not go to —00 in reality. However, in amplifiers, the system is often set up such that some of the modes of the 
system are just this type of negative energy modes at least for small enough perturbations of the system. In a laser for 
example, pumping the atoms to their excited state gives a systems whose small fluctuations are exactly the above sort 
of negative energy harmonic oscillator. At large amplitudes, those modes will saturate and the system will become 
non-linear (when a significant portion of the atoms have made the transition from the excited state to the ground 
state). Thus this model is not a good model for all regimes of such an amplifier, but is a good approximation as long 
as one is concerned only with its small signal behaviour. Note that that small signal regime can be one in which the 
excitations of the ip field are much much larger than the size of the mean quantum or thermal noise in that field. 




In(Freq) 

FIG. 1: Figure 1. Amplification for the Q — situation. Note the similarity to the amplification of a standard integrated 
circuit amplifier. 



The annihilation operators for a specific mode 4>i(t — x), assumed in the distant future to be far from the origin 
x = is 



.4) 



b*U^(t,x) - d t <f>i(t - x)*$(t,x)) dx 



(56) 



where $ and 11$ are the quantum field and conjugate momentum operators in the Heisenberg representation. The 



which, if (hi is normal- 



annihilation and creation operators obey the usual commutation relation [a^al.] = i(<f>i, 

ized, is the usual commutation relation for annihilation operators. Similarly, the annihilation operator for a tp mode 



iPj 



Hi 



• J j • 



$) =i 



j (V>JII*(i, x) + d t ipj(t - *)**(*, a;)) dx 



The relation between the outgoing modes of the ip and (j> fields to the ingoing are then 



-u? + n 2 - im{e 2 + e 2 ) 
-uj 2 + fl 2 - iuj(e 2 - e 2 ) 



' a i>,ui,in^ 1 



+ n 2 -iw{e 2 -i 2 ) 



= A u 



where the amplification factors obey 



B u a 



^^■tp,u},in 



\A U 



\BJ' = 1 



(57) 

(58) 
(59) 

(60) 



as required. 

(2 ~2 \ ^ 
e 2 _~ 2 ) and thus to get a large amplification 

we require that e 2 — I 2 « e 2 + e 2 . 

An interesting situation occurs if we take the central oscillator to be a free particle, so that $7 = 0. The amplification 
is then roughly constant for a frequency range around lo = until uj rj e 2 — e 2 and then falls at 6dB/octave for 
frequencies higher than that until ui ~ e 2 + e 2 at which frequency the amplification goes to 1. This is just the 
behaviour one has for many amplifiers- eg the gain curve of a transistor amplifier. 

The quantum noise temperature for the O = case (assuming both of the input modes are at zero temperature) 
and assuming e 2 — e 2 << 2e 2 is given by 



T = — 



In 

' 4e 



( \B\ 2 



111 



2 + (e 2 +e 2 ) 2 



In 



;w << 2e 2 
;w » 2e 2 



(61) 



(62) 



Ie, the quantum noise temperature diverges near uj = 0, and also for large uj and achieves a minimum at around 
uj ~ 4e 2 . For uj smaller than that this value, the temperature is larger than uj, but for larger uj, the temperature is 



smaller than u> and the number of particles in those modes will be much less than unity (ie, one is in the Wien tail 
of the distribution) . 

Note that I have assumed that both the output and input channels are the <fi field. One could also choose the input 
channel to be the <$> field and the output the ip field in which case B^ would be the amplification factor instead, and 
the amplification would fall to zero, rather than 1 for very large frequencies. 

What is clear is that this simple model for a continuum amplifier captures many of the features of a real amplifier 
and can be applied to a wide variety of amplifiers [7[. 

III. BLACK HOLES 

One of the more fascinating forms of amplifier is that provided by a black hole[8|. Hawking 9] showed that that the 
relationship between the ingoing and outgoing modes of black hole could be written as 

e /3u/2 e -/3w/2 

4>u;,out = -j====(f> uin + — j=====(f> uin (63) 

y/smh(puj/2) ^/sinh(paj/2) 

where ipuj.out is a positive norm mode escaping from the black hole with frequency u>, and 4^ in * s an m °de ingoing 
toward the star which will eventually form the black hole made up entirely of positive norm positive frequency modes 
of the ingoing the field, and cj>~ in are modes made up entirely of the usual negative norm, negative frequency ingoing 
modes. The relation between the ingoing and outgoing modes is unusual in that the relation between the ingoing 
energies and outgoing energies is bizarre. 

Let us take as a model for black hole formation the collapse of a spherically symmetric null shell of dust. Before 
the collapse, spacetime is flat, with null coordinates U and V and with "radius" (-^t of the circumference of the 
spheres of spherical symmetry) given by r = (V — U)/2. Let me choose the coordinate v so that V — corresponds 
to the shell of dust. Outside the shell the metric is Schwartzschild, with null coordinates u, v with (v — u)/2 — 
(r — 2M) + 1M ln C^M^ )i an d with the null shell given by v = 0. Along the null shell the requirement that the 
circumferential radius be continuous across the shell gives us the relation between the U and u coordinates as 

u = (+£/ + AM) - AM In ( _ ( u + * M ) \ ( 64 ) 

V 4rlW J 

Thus, if we have a wave-packet of the form of 

<f>out = S{u)e iua (65) 

where S is a relatively slowly varying envelope concentrated around u » AM we will have that the form for the 
incoming wave function will be 

4MV) » S (4MM^2) (^)~ l4MU ;V<~4M (66) 

;V>-AM (67) 

This can be written in terms of the ingoing positive norm modes which are pure linear combinations of the ingoing 
positive norm modes e~ lUV with ft > 0. 

{^WY^ ;V>-4M (68) 



-\/sinh(47rAfw) 

e - 2 "M» f-iM-V\ i4MuJ 



; V < -AM 



%,u 



^ { 4M±vy^ v> _ m ^) 

where in each case uj > 0. These two positive norm modes correspond to the a, b incoming modes discussed in the 
first section. 

Note that in this case the two types of mode a, b correspond to different types of the same incoming field <f>i n . The 
frequency of the ingoing mode which goes as 4 | f -,_ i4M ^ is approximately 

o i i ({4M + VY HMu \ 4Mo; u/iM 

\l fa idvln ~ ~ we ' (70) 

V \ \ AM I AM + V V ' 
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Ie, the frequency of the incoming mode which creates an outgoing mode of frequency ui at retarded time u is exponential 
in that retarded time. For example for a retarded time 1 second after a solar mass black hole forms, the incoming 
frequency corresponding to an outgoing frequency of u is about cue 10 which is e 10 times a frequency corresponding 
to the mass of the whole universe. Thus the energy of the incoming modes which are amplified with an amplification 

2Moj — u 

factor of e by the black hole amplifier, have their energy decreased by a factor of eiw. Amplification 

\J2 sinh(47rMaj) 

does not imply energy amplification. The black hole, as an amplifier, amplifies the amplitudes (norms) by a thermal 
amplification factor, but de-amplifies the energy by a term with is exponential in the time after the black hole forms. 
It is this feature of the black hole as an amplifier that makes it unique. 

IV. DUMB HOLES 

In 1981 lpi| suggested that many of the features of the black hole particle creation could also be captured by 
what I have since called Dumb holes- analogs in condensed matter system which have horizons and mimic many of 
the features of black holes, including the output of quantum noise as the analog of Hawking radiation. In the case 
of Hawking radiation, if one traces back the the emitted radiation into the past, it is squeezed against the horizon 
exponentially until one gets back to the time when the black hole was originally formed. Only then can the backward 
propagating modes escape toward infinity, but with absurdly high frequencies and short wavelengths. 

Dumb holes (named after the original usage of the term to mean "unable to speak" , and not the more modern usage 
to mean "stupid"- ie, "Stumm" not "Dumm" in German, or "Muet" not "Stupide" in French), are systems in which 
the flow of some background material is such that some wave in that material, which has a non-trivial dispersion 
relation, to experience a horizon, a surface through which the low frequency components of the wave cannot travel 
from inside to out because the ingoing velocity of the material equals the outgoing velocity of the wave. Such systems 
experience quantum noise in the same was as black holes do- ie they emit a thermal spectrum of radiation of that 
quantized wave whose temperature is determined, instead of by the mass as in the black hole case, by the nature of 
that background flow at the horizon. 

In the case of dumb holes, the same thing happens as for black holes initially (finally since we are tracing the 
modes backward in time?) and the backward-in-time modes are exponentially squeezed against the horizon. But at 
sufficiently short wavelength the dispersion relation and in particular the group velocity of the modes changes and 
the mode escapes from the horizon with very short wavelengths. Depending on the nature of the dispersion relation 
the outgoing mode can either be dragged in from large distances because the group velocity is now much less than 
the velocity of the fluid far from the horizon in the case where the dispersion relation makes the group velocity small 
at high frequencies, or can now travel out from inside in the case in which the dispersion relation has a group velocity 
at high frequencies much larger than the velocity of the fluid inside the horizon. 

If we assume the flow to be stationary, the waves see a time independent situation, and, while the wavelength can 
change drastically, the frequency is conserved in the lab frame ( but not in the fluid frame). In Figure 2 we have 
a graph of such a dispersion relation in a still fluid, in which the group velocity falls at high frequency. If the fluid 
moves with some velocity which is smaller than the long wavelength speed of the wave in still water, the dispersion 
relation looks as in figure 3, while if the fluid is moving with a velocity higher than the long wavelength speed of the 
wave, figure 4 gives the dispersion relation. 

As an example, let us assume we have a 1+1 dimensional wave in a fluid whose still fluid dispersion relation is 

u? = F 2 {k) (71) 

In the moving fluid the dispersion relation will be 

w = vk ± F(k) (72) 

where the two branches represent the left and right moving waves. We will be interested in the minus sign which will 
represent waves which are trying to move against the flow. 
The Lagrangian for such a fluid could be given by 

£= f{dt(t>-v{x)d x (t>Y-{F{id x )(t>Y (73) 



where I have assumed that F is an odd function of k. 
The norm is again 



'/ 



% / 4>*{d t 4> - v{x)d x 4>)dx + CC (74) 




FIG. 2: Figure 2. Dispersion relation for a "subluminal" case in still fluid. (This is in fact the relation for surface waves on 
water, with to — ^J gk tanh(fch)) 




FIG. 3: Figure 3. Dispersion relation for waves as in figure 1 with the flow rate being 'subluminal" (ie, slower than the velocity 
of the waves at zero wave- vector) . The horizontal line is for frequency 1, and the three possible wave- vectors are kg which is a 
positive norm wave, both of whose phase and group velocity is to the right, kf has positive phase velocity and negative group 
velocity, and has positive norm, k^ has negative phase and group velocity and negative norm 



If we assume v to be a constant, and <j)(t, x) — fy^e then the norm will be 

< <p, <p >oc 2(w- vk)\(j> k \ 2 



(75) 



The important point to note is that the norm is positive or negative, depending on the sign of u> — vk Since u> is given 
by the above dispersion relation, we have 



<0,0>=2OF(fc))W 



(76) 
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FIG. 4: Figure 4. Dispersion relation for waves as in figure 1 with the fluid flowing faster than the velocity of the still water 
waves at small wave-number. There is only one possible wave-vector ki for the given frequency, and it has negative phase and 
group velocities and negative norm. 



Ie, the sign of the norm depends not on the value of ui but the sign of the still water dispersion function for that 
mode. Modes with positive u> can have negative norm, and modes with negative u> can have positive norm. 

Let us imagine that we have a flow where the fast flow occurs to the right and the slow flow to the left. There will 
be a horizon between the two. Now consider modes with the frequency indicated in the diagrams. The modes which 
have group velocity toward the horizon are the modes with wave- vectors k\ and hi. We see, that just as in the above 
continuum model of an amplifier, the ingoing modes with a given positive frequency u> come in two flavours, the ones 
with positive norm (ho and </>o) and ones with negative norms, but the same positive frequencies (k^ and ipo). The 
outgoing modes- travelling away from the horizon, or away from the oscillator in the continuum model- also come in 
the same two flavours, the positive norm modes (fci and <j} ou t) and the negative norm outgoing modes (&2 and tpout)- 
In the continuum model the coupling between these modes which leads to amplification is the harmonic oscillator. In 
this dumb hole model, the coupling between the modes is provided by the non-adiabatic, spatially dependent changes 
in the background flow given by the changing velocity v{x). 

In the case of the dumb hole, but not in the oscillator coupling, the effective temperature of the emitted quan- 
tum noise from this amplifier has a constant temperature, independent of frequency, at least for low frequencies. 
This differs significantly from the continuum amplifier mentioned above where at low frequencies the temperature 
diverges. [13 [Il|[i3 



APPENDIX 



To show that the C state is left in a thermal density matrix after tracing out over the D states, it is easiest to do 
so if we assume that the state of the system is the vacuum state for the a and b inputs. In this case the condition 



«|0) ob = b\0) ab 



(77) 



becomes 



which has solution 



cosh(60C + sinh(>).D t |)) ab =0 
cosh(6»)D + sinh(/x)C t |)) a6 - 



(78) 
(79) 



|0)„fc 



c ,tanh(/i)C t .Dt 



|o> 



CD 



(80) 
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where |0)c_d is the state annihilated by the C, D operators. Tracing over D by using the quanta eigenstates \m) — 
( , ' , 10) d we have 

Vm! ' ' 

Tr D \0} ab (0\ ab = J2(m\ D ^ tanh( M ) r ^ ^ Y \0}cd(0\cd J2 tanh( M )^^^|m) g (81) 

m r s 

= ^tanh( Ai ) 2m i 7 i r |0)c(0|c^= (82) 

ym\ \/m\ 

= ^ e m(21n(tanh( ' i)) |m>c(m|c (83) 

771 

_ e 21n(tanh( AI ))C t C fg^s 

which is a thermal density matrix with thermal factor A = — 21n(tanh(/i)).|8j 
It is also clear that if one began with the two mode squeezed state 

|^ =e -21„(ta„h( AI ))at b t| ^ & (85) 

That state in terms of the output modes would just be the vacuum state |0)cd which would minimize the output 
noise. 

To show that if the input state is a thermal state in each of the channels, 

p = e A ^ a e h ^ b (86) 

the output state of the C channel is also a thermal state, I found it easiest to go use path integrals. Using 

flt + a .a^ + a 

*=- ?r fl.—^ (87) 

&t + b tf + b 

we can write the density matrix, 

p = ^e-*Ao( P =+x 2 ) ( 89 ) 

between the initial and final eigenstates of x,y operators as We can take as our example the density matrix p = 
e s(Px+ x to get 

f 1 1 

W)p{x\ = / {x,y\\p xN ){p xN (l - -A a (p2 + x 2 )/N\x N _ 1 ) ( Pxi \(l - -A a (pl + x 2 )/7V|^_i)... 

....(p xl \-A a {pl + x. 2 )/N\x)tt k dp xk dx k (90) 

w e T.^(^^)-^l^)^ nkdpxkdXk (91) 

= feIo lp ^ T)±{T) -^ MT)2+x(T}2)dT U T dp x (T)dx(T) (92) 

Completing the squares in the exponent with respect to p x and doing the p x integrals we get 

- 1 f 1 (± 2 
(x'\p\x) = J e-^ — + Ax 2 )dT\U T dx{r) (93) 

where x(l) = x' and x(0) = x. 

Similarly for the two modes, we have 

(x',y'\p\x,y)=N [\ e i S ¥+ A ^M 2 +^+ A «(^ 2 )^n r fe(T)^(r) (94) 



ii 



where the path integral is taken over all paths x(r),y(r) such that x(0) = x, y(0) — y, x(l) = x', y(l) = y' ■ 
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As usual we can do a change of variables of the path integral, such that 

x(r) = x(t) - X(t) (95) 

y(r) = y(r) - Y(t) (96) 

where X(t), Y(t) obey 

X = A 2 a X (97) 

Y = A 2 b Y (98) 

and where X(0) — x, Y(0) — y,X(l) = x', Y(l) = y'. The boundary condition on the tilde variables is 

x(0) = x(l) = j/(0) = j/(l) = (99) 

The exponent of the path integral then becomes 

s = \ J! ^a7 + AaX{T)2 + ^aT + A ^ T ^ dT ( 10 °) 

1 f X{T)2 + K a X(rf + 1^1 + K b Y{rf) + \ f *££ + K a x{rf + M- + K b y(rf) (101) 



2 J A a w A 6 " w "2j A a w A„ 

where the cross terms between X, 5 and Y, y vanish by integration by parts and because X, Y obey the equations of 
motion, and because x, y are zero at the endpoints. 

Doing an integration by parts on the X, Y terms, and using that they obey the equations of motion we get that 
the only contribution to the integral is from the endpoints where we get for the integrand 

S = -\{{x'X{l) - xY(0))/A a + (y'Y(l) - yY(0))/A 6 )) (102) 

~\j^K7 + KHt)2 + ^at + A »y^ dT ( 103 ) 

where the path integral now is over paths where the endpoints of the tilde variables are all 0. The contribution of 
the second part (the integration over the tilde variables) to the path integral is independent of the values at the end 
points x, x' ', y, y' so it simply multiplies the path integral by a constant which can be absorbed into the normalisation 
factor N. The solution for X, Y with the given boundary conditions is 

x = / inh(A(l-r)) + ^sinh(Ar) 

sinh(A a ) sinh(A ) 

y = / inh(A,(l-r )) + sinl^r) 

smh(A{,) sinh(Afe) 

which gives as the only non-trivial contribution to the integrand 

p{x'y';xy) ex e s (106) 



with 



~ = 1 x'X(l)-xx(0) | y'Y(l)-yY{0) 



2 A„ A b 

sinh(A a ) ' ' ' " ''" sinh(Ab) 



-- ( (x 2 + x' 2 ) coth(A a ) - 2x3;' - — - + {y 2 + y' 2 ) coth(A b ) - 2yy' . ) (108) 



and 



p(x'y'; x, y) = (x'y'\p\xy) = Ne s (109) 

Now, we want to take the trace of this density matrix over all D states. Defining 

z-£+£. ,iio, 
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we have that 

\x, y) = | cosh(/i)Z - smh{p)W, cosh{p)W - sm(p)Z) (112) 

and the trace of p over D becomes 
(Z'\ Tr DP \Z) = 

p{cosh(p)Z' - sinh(yu)W, cosh(/i)VF - sm(p)Z'; cosh(p)Z - sinh(p)W, cosh(^)VF - sm(p)Z)d\fL13) 



Since the integrand is an Gaussian exponential in the three variables Z, Z', W, after the integration over W, (since 
the coefficient of W 2 is independent of Z, Z') the result is also Gaussian in Z and Z' and is symmetric in Z, Z' . Ic, 
it is also a thermal state. Explicit calculation, by completing the squares in the exponent of the integrand for W, 
shows it is of the form 



(Z'\Tr D p\Z) ex e ~ ( "°MA C )( z + z )ssm^ zz ) 



(114) 



which is again a thermal density matrix with thermal factor Ac. While one could actually evaluate the terms in order 
to determine what Ac is in terms of p, A a , A;,, it is far easier to do this by the procedure in the main section and 
simply evaluate Tr(pC Ji C) to determine the thermal factor. 
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